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ABSTRACT 

We present in this paper the PolEMICA (Polarized Expectation-Maximization Independent 
Component Analysis) algorithm which is an extension to polarization of the SMICA (Spectral 
Matching Independent Component Analy sis) temperature multi-de tectors multi-components 
(MD-MC) component separation method dDelabrouille et al. 20031) . This algorithm allows us 
to estimate blindly in harmonic space multiple physical components from multi -detectors po- 
larized sky maps. Assuming a linear noisy mixture of components we are able to reconstruct 
jointly the anisotropics electromagnetic spectra of the components for each mode T, E and B, 
as well as the temperature and polarization spatial power spectra, TT, EE, BB, TE, TB and 
EB for each of the physical components and for the noise on each of the detectors. PolEM- 
ICA is specially developed to estimate the CMB temperature and polarization power spectra 
from sky observations including both CMB and foreground emissions. This has been tested 
intensively using as a first approach full sky simulations of the Planck satellite polarized chan- 
nels for a 14-months nominal mission assuming a simplified linear sky model including CMB, 
and optionally Galactic synchrotron emission and a Gaussian dust emission. Finally, we have 
applied our algorithm to more realistic Planck full sky simulations, including synchrotron, 
realistic dust and free-free emissions. 



Key words: 

analysis 



Cosmic microwave background - Cosmology: observations - Methods: data 



X 



1 INTRODUCTION 

Mapping the Cosmic Microwave Background (CMB) polarization 
is one of the major challenges of future missions of observational 
cosmology. CMB polarization is linear and therefore can be de- 
scribed by the first three Stokes parameters I, Q and U which are 
enerally combined to prod uce three fields (modes), T, E and B 
Zaldarriaga & Seliak 19971) . The polarization of the CMB photons 
carries extra physical informations that are not accessible by the 
study of the temperature anisotropies. Therefore its measurement 
helps breaking down the degeneracies on cosmological parame- 
ters as encounter with te mperature anisotropies measurements only 
dZaldarriaga et al. 19971 ). Furthermore, the study of the CMB po- 
larization is also a fundamental tool to estimate the energy scale of 
inflation which has been proposed to solve the problems of flatness, 
of isotropy and of the seed perturbations for the formation of the 
structures in the Universe. Inflationary models predict the presence 
of tensor perturbations of the metric which will lead to an unique 
signature in the CMB polarization B modes. The detection of the 
latter would be a strong proof of such an epoch and also a way to 
constrain the energy scale at which inflation occu rs by measuring 
the tensor to scalar ratio, r ^Turner & White 1996b . 

Since the beginning of the CMB anisotropies ob- 
servations with the COsmic Background Explorer 
(COBE) iSmoot et al. 19921) . a great amount of experi- 
ments have been designed to determine the CMB tem- 
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contrast, the polarization anisotropies, which are between 2 
and 5 orders of magnitude weaker than temperature ones are not 
accurately measured yet. A first detection of th e CMB E modes ha s 
been perfor med bv DASI Jkovac et al. 20021; iLeitch et al. 2005]) , 
CAPMAP jBarkats et al. 2005b . CB I (IReadhead et al. 20041) 
and mor e recently by B OOMERanG l lMontrov et al. 20051) and 
WMAP JPage et al. 2006b . The TE temperat ure-polarization cross 
correlation has been measured by W MAP JPage et al. 2006b and 
BOOMERanG JPiacentini et al. 2005b . No detection of the CMB 
B modes has been reported yet. Nevertheless, constraints on the 
tensor to scalar ratio, r, have been set by the WMAP team. T hey 
set an upper limit of r < 0.55 (95% CL) dSpergel et al. 2006b for 
the temperature a nd polarization analysis and of r < 2.2 (95% CL) 
dPage et al. 2006b for a polarization-only analysis. 

The detection of such low signals is possible by improving 
the instrumental sensitivity, but this is not the only issue in the 
determination of the CMB polarization power spectra. Other as- 
trophysical emissions as for example the diffuse Galactic emis- 
sion including free-free, dust and synchrotron and the extragalactic- 
sources emissions also contribute to the sky brightness at the fre- 
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quencies of interest for CMB studies, and therefore must be ef- 
ficiently subtracted. These foregrounds are particularly important 
for the study of the CMB polarization. Excluding the free-free 
emission which is not polarized, the other contributions are ex- 
pected to be significantly polarized with similar power on the E 
and B modes. Recent measurements of the Galactic sync hrotron 
polarization emissio n at 1.41 GHz i Wolleben et al. 2005b and at 
23 GHz by WMAP dPage et al. 2006h~ show this emission is sig- 
nificantly polarized at large angular scales. Further, Archeops 
measurements at 353 GHz show that the Galactic dust diffuse 
emission is pol arized up to a lev el of 5 to 10 % both in the 
Galactic center l|Benoit et aL 2004h and at high Galactic latitudes 
dPonthieu et al. 20051) . Finally, for the polarization of extragalactic 
point sources the sparsity of the data available make s reliable pre- 
dictions difficult dTucci et al. 20041 ; iHildebrand 19961) . 

A direct subtraction of these foreground contributions on 
the CMB data will require an accurate knowledge of their spa- 
tial distributions and of the electromagnetic spectra of their 
anisotropies. For the synchrotron emis sion a full sky m ap at 
408 MHz in temperature is available dHaslam et al. 19821) and 
more recently the WMAP team provided a map at 30 GHz 
from the MEM dec omposition of the first year observations 
jBennett et al. 2003bT) . A fak e polarized synchrotro n emission tem- 
plate was construct ed by dGiardino et al. 20021) based on the 
Parkes 2400 MHz dDuncan et al. 19971) and Haslam 408 MHz 
dHaslam et al. 1982b surveys. Furthermore, the electromagnetic 
spectrum of synchrotron anisotropies and its spatial distribution 
are neither accurately known in temperatur e nor polarization al - 
though a first estimate was produced by dGiardino et al. 20021) . 
Recently, the 23 GHz polarized WMAP data is u sed as a tracer 
of the synchrotron polarization dPage et al. 20061) . For the ther- 
mal dust emission a full sky map at 100 fim as well as tem- 
plates for CMB use were ex t racted from the IRA S and FI- 
RAS data dSchlegel et al. 19981: iFinkbeiner et al. 19991) . No real- 
istic template exists for the dust polarized emission, although a 
fake one, based on the polarization angles measured at 23 GHz 
has been constructed by the Planck collaboratiorQ The dust emis- 
sion in temperature can be approximated by a grey body spectrum 
of mean temperature 1 7 K and emissivity between 1.7 and 2.2 
dFinkbeiner et al. 19991 : lLagache et al. 20031) . Currently no mea- 
surement on the electromagnetic spectrum of the dust polarized 
emission is avai lable although it is expected to be the same that 
for temperature l ljones et al. 19921) . 

To try to overcome the above limitations, a great amount 
of work has been dedicated to design and implement algo- 
rithms for component separation which can discriminate be- 
tween CMB and foregrounds. These methods can also ex- 
tract, directly from the CMB data, the emission properties of 
foregrounds. Wiener filtering has been successfully tested as- 
suming known Gaussian priors for each component and with 
the electromagnetic spectrum of the anisotrop ies as an input 
dTegmark & Efstathiou 19961: iBouchet et al. 19991) . Maximum en- 
tropy based methods (MEM), assuming entropic priors for the spa- 
tial distribution of each of t he component, have been intensively 
used for small sky patches dHobson et al. 19981) and extended to 
full sky analysis dStolvarov et al. 2002b . They were adapted to 
account for spatial a nisotropies in t he electromagnetic spectra 
dStolvarov et al. 20 04). More recently, dEriksen et al. 20051) has de- 
veloped a new method to perform CMB component separation by 
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parameter estimation and applied it to temperature simulations of 
the Planck satellite experiments. Independent Component Anal- 
ysis (ICA) techniqu es have also been applied to Planck simula- 
tions inj£mp_eraturej Maino et al. 20021) and ex tended to polariza- 
tion dBaccigalupi et al. 2004l ; IStivoli et al. 20061) using the FastICA 
algorithm. These methods require no prior on the spectral or spatial 
distribution of the components but can not make use of the avail- 
able physical knowledge on the foreground and CMB emissions. 
In addition, th e Spectral Matching Inde pendent Component Analy- 
sis (SMICA) dDelabrouille et air20 03) has been developed to con- 
sider both the fully blind analysis for which no prior is assumed 
and the semi-blind analysis incorporating previous physical knowl- 
edge on the astrophysical components. This algorithm, based on th e 
Expectation-Maximization algorithm (EM) dDempster et al. 1977b . 
uses the spectral diversity of the components and was developed for 
temperature only. We present in this paper, PolEMICA (Polarized 
Expectation-Maximization Independent Component Analysis), an 
extension of this method to polarization including both the blind 
and semi-blind analysis. 

This paper is organized as follows. A simple model of the 
microwave sky emission in temperature and polarization is de- 
scribed in section [2] Section [3] presents the PolEMICA multi- 
detectors multi-components (MD-MC) blind component separation 
algorithm. Section|4]describe the simulations of the Planck satellite 
experiment used for testing the algorithm. We present in section [5] 
the application of PolEMICA to the Planck simulations with a sim- 
plified model to test the algorithm's performances. Finally, in sec- 
tion[6] we apply our algorithm to more realistic Planck simulations 
and discuss the separability problem in this case. We summarize 
and conclude in section]?] 



2 MODEL OF THE MICROWAVE AND SUB-MM SKY 

2.1 Multi-detectors Multi-components model 

To constrain cosmological models, CMB experiments have to reach 
an accuracy which is well below the expected level of contamina- 
tion from astrophysical foregrounds, in temperature and even more 
critically in polarization. Therefore, an efficient separation between 
CMB and foregrounds is crucial for the success of future polariza- 
tion experiments. To perform such a separation, the diversity of 
the electromagnetic spectra of the anisotropies and of the spatial 
spectra of the components is generally used. Observations from a 
multi-band instrument can be modeled as a linear combination of 
multiple physical components leading to what is called a Multi- 
Detectors Multi-Components (MD-MC) modeling. 

Assuming an experiment with n y detector-bands at frequen- 
cies y, and n c physical components in the data, for each Stokes 
parameter (/, Q and U) and for each pixel on the sky map we can 
write 



/(r) = A v c ■ s c (r) + n v (r) 



(1) 



where s c (r) is the map of the c component, zi v (r) refers to the noise 
map for each v, band and A v c which is called the mixing matrix, 
gives the electromagnetic spectrum behavior for the component c 
and frequency v, . Beam smoothing and filtering effects are not con- 
sidered in this work. 

As in the temperature case, it is more convenient to work in 
spherical harmonics space, where equation (QJ can be rewritten 
independently for each [l, m] assuming a full sky coverage. Thus 
equation (1} reads for X = [T,E,B] and for each frequency band 
and for each {€, m) 
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V* /li'.X c,X , v,X / os 

= Zj A '' s (>» (2) 

c=l 

where yzf is a vector of size (3 • n v ■ ni ■ n m ), s c e "* is a (3 • n c ■ ni • n m ) 
vector and n v : is a vector of the same size than j£ . A v c ' x is a matrix 
of (3 • n v ) x (3 • n c ) elements formed from the mixing matrix of each 
of the modes, T, E and B. 

The aim of the component separation algorithm presented in 
this paper is to extract A* , s c ( ' m and nj£ from the yj£ sky observa- 
tions. 



v(GHz) 


30 


40 


70 


100 


143 


217 


353 


CMB 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


Sync. 


1.0 


0.46 


0.11 


0.045 


0.021 


0.012 


0.014 


Dusl 


0.0006 


0.001 


0.003 


0.008 


0.021 


0.088 


1.0 


Free-free 


1.0 


0.56 


0.19 


0.10 


0.061 


0.046 


0.071 


Noise / 


4.12 


4.03 


4.06 


1.47 


1.0 


1.47 


4.54 


Noise P 


2.91 


2.95 


2.98 


1.51 


1.0 


1.51 


4.59 



Table 1. Electromagnetic spectrum of the anisotropies, in arbitrary units, of 
the CMB, dust, synchrotron and free-free components at each of the Planck 
channels. Relative noise levels in the temperature and in the polarization 
maps taking as reference the 143 GHz channel are also displayed. 



2.2 Simulated microwave and sub-mm sky 

Following the MD-MC model discussed above and given an ob- 
servatio nal setup, we const ruct, using the HEALPix pixelization 
scheme dGorski et al. 19991) and in CMB temperature units, fake 
/, Q and U maps of the sky at each of the instrumental frequency 
bands. For these maps we consider three main physical components 
in the sky emission: CMB, thermal dust and synchrotron. Concern- 
ing emissions that are not expected to be significantly polarized, 
we have optionally considered unpolarized free-free emission and 
not taken into account the SZ emission as we are interested in dif- 
fuse emissions. When no free-free is considered and in any cases 
for the SZ emission, we suppose that they have been successfully 
removed from our sky maps in temperature. We have assumed that 
the free-free emission is not polarized. Instrumental noise is mod- 
eled as white noise. 



CMB 

The CMB component map is randomly generated from the polar- 
ized CMB angular pow er spectra computed with the CAMB soft- 
ware dLewis et al. 2000h for a set of given cosmological parame- 
ters. In the following we have used Ho = 71 km ■ s" 1 • Mpc~ , Cli, = 
0.044, Q m = 0.27, n A = 0.73 and t = 0.17 that are the values of 
the cosmolog ical concordance m odel according to the WMAP one 
year results dSpergel et al. 2003). We also consider gravitational 
lensing effects as described in Jchallinor & Lewis 200alHu 2000l : 
IChallinor & Chon 20021 ; Ibkamoto & Hu 2003h . 



Synchrotron 

For the diffuse Galactic synchrotron emission we use the 
tem plate maps in temp erature and in polarization provided 
by dGiardino et al. 20020 . These template maps were derived 
in temperature dir ectly from the Haslam map at 408 MHz 
faaslam et al. 1982b . The polarization maps in Q and U were 
constructed from the intensity map from a constrained real- 
ization of the polarization angles using the Parkes 2400 MHz 
survey jDuncan et al. 1997h . A template of the spatial varia- 
tions of the synchr otron spectral index is also provided by 
dGiardino et al. 20021) . Here we have chosen to use a constant spec- 
tral index equal to the mean of the spectral index map, a = -2.77, 
so that the simple linear model of the data holds. A more realistic 
treatment of the synchrotron emission will require a specific model 
of the spectral index spatial variations in order to ensure separabil- 
ity between components. 



Dust 



★ Simplified-d«.?f model: We have first considered a simpli- 
fied dust mo del which is Gauss ian and derived from the power-law 
model from dPrunet et al. 1998!) to describe the dust angular power 
spectra in temperature and in polarization. This model, although not 
fully realistic, is not spatially correlated to the synchrotron emis- 
sion and helps us to extensively test the properties of the separation 
method. We have renormalized this model to mimic at large an- 
gular scales the TE cross power s pectrum measured by Archeops 
at 353 GHz dPonthieu et al. 20051) . The rms of the final dust map 
is probably overestimated as we do not account for the variation 
of the dust emission with respect to Galactic latitude. The power 
spectra models are computed at 100 GHz in /jK rj units. /, Q and 
U full-sky maps are generated randomly from these power spectra. 
We extrapolate them to each of the frequency of interest by assum- 
ing a grey body spectrum with an emissivity of 2. Finally, the yuKm 
maps are converted into yuK C MB units. 

★ Realistic-diw/ model: Seco ndly we have used the Planck Sky 
Model 1 polarized dust template i Baccigalupi 20031). It was mod- 
eled using model number 7 of dFinkbeiner et al. 19991) . This model 
is normalized to the IRAS 100 fim emission map produced by 
dSchlegel et al. 1998h . For polarization, a constant polarization de- 
gree of 5% is assumed, and the same polarization angles than for 
the synchrotron model are used. For both temperature and polariza- 
tion, we assume a grey body emission with an emissivity of 2. 



Free - free 

The free-free component is derive d from the Planck Sky Model 1 . It 
is based on the H-a template by dDickinson et al. 20031) . We have 
assumed a constant spectral index of -2.1 in Rayleigh- Jeans units. 
Expected to be no signi ficantly polarized ex cept in particular HII 
regions (less than 10% dKeating et al. 19981) ). the free-free emis- 
sion <2 and U maps are set to zero. 



Noise 

Noise maps for each channel are generated from white noise 
realizations normalized to the nominal level of instrumental noise 
for that channel. 

The electromagnetic spectra of the anisotropies, in arbitrary 
units, for the CMB, dust and synchrotron emissions are displayed 
in table[T](these are the values used in the mixing matrix A) for the 
Planck satellite simulations presented in section|4] We also present 
the relative noise level taking as reference the 143 GHz channel. 
The noise levels used at 143 GHz are 6.3 /jKcmb (in temperature) 



We have used through this article two different dust models. 



1 [httpTT/www . cesr.fr/~bernard/PSM/ 
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and 12.3 /uK C mb (in polarization) per square pixels of side 7 arc min 
and for a 14-months Planck mission (Pla nck Consortium 20051) . 



3 A MD-MC COMPONENT SEPARATION METHOD 
FOR POLARIZATION 

3.1 MD-MC model for the temperature and polarization 
power spectra. 

To reduce the number of unknown parameters in the model de- 
scribed by equation l[2j, it is interesting to rewrite this equation 
in terms of the temperature and polarization auto and cross power 
spectra. This will considerably reduce the computing time with no 
loss of information. 

We define the density matrices associated with the data, y, the 
physical components, s and the noise, n, as follows 



2( 



. m 

— V v'- x 
+ i Zj Af " 



X i'Xl 
Xf m 



(3) 



where i represents frequency, v, for the data and noise matrices 
and component, c, for the physical-components matrix. Averaging 
over bins on t we obtain 



"b o-rr\ 



where Dt is the set of t values which contributes to bin b and nt, is 
the number of such multipoles. In the following, N b , represents the 
total number of bins used in the analysis. 
For each bin b, equation $2^ reads 



R,.(b) = AR s (b)A' + R„(b) 



(5) 



where R y (b) and R„(b) are (n Y ■ 3) x (n v ■ 3) matrices and R s (b) is a 
(n c ■ 3) x (n c ■ 3) matrix. 

To fully understand the component separation algorithm de- 
scribed below it is interesting to have a closer look to the content of 
the three density matrices defined above (see appendix for a con- 
crete example). 

R y (b) represents the input density matrix computed from the 
observed multi-band data. This matrix is composed of 3 x 3 sym- 
metric blocks each of them containing in the diagonal the auto- 
power spectra, TT, EE, BB and in the off-diagonal the cross-power 
spectra TE, EB and TB. A single block represents either the auto- 
correlation of a single channel (for diagonal blocks) or the cross- 
correlation between two channels (for off-diagonal blocks). 

Assuming that the physical components in the data are statis- 
tically independent and uncorrelated makes the R s (b) a 3 • n c x 3 • n c 
block diagonal matrix. As above each block, corresponding to the 
f lh physical component, contains in the diagonal the auto-power 
spectra, TT, EE, BB and in the off-diagonal the cross-power spec- 
tra TE, EB and TB. 

We also assume that the noise is uncorrelated between chan- 
nels and therefore, R„(b) is a diagonal matrix containing the noise 
auto power spectra TT, EE and BB for each of the channels. 



3.2 Spectral matching algorithm 

From equation l[5) we observe that, for each bin b, the data density 
matrix R y (b), of size 9 x n\, is fully defined by the set of parameters 
9(b) = [A,R s (b), dmg(R„(b))} which corresponds to a total of 



3 x n v x n c + 6 x n c + 3 x parameters. This indicates that from 
the CMB data set and under the hypothesis presented above it is 
possible to simultaneously estimate the mixing matrix, the physical 
component's temperature and polarization power spectra and the 
noise's temperature and polarization power spectra for each of the 
channels. Further, assuming white noise in the maps only three 
noise parameters per channel (for TT, EE and BB) need to be 
estimated for the entire range in £. This reduces the overall set of 
parameters to 3 x n v x n c + 6 x x N b + 3 x n v where N b is the 
number of bins. 



The I ikelihood function 

To estimate the above parameters from the data we have extended 
to the ca se of polarized data the sp ectral matching algorithm devel- 
oped by dDelabrouille et air2003Q for temperature only. The key is- 
sue of this method is to estimate these parameters, or some of them 
(for a semi-blind analysis), by finding the best match between the 
model density matrix, R y (b), computed for 6(b) and the data den- 
sity matrix R y (b) obtained from the multi-channel data. Assuming 
that the different physical components and the noise are realiza- 
tions of Gaussian stationary fields (Wittle approximation), the log- 
likelihood function of the form 



Nb r 

(4) -X = J] n >- Tr{Ry(b)R- y \by} + log det 



(6) 



is a reasonable measure of the mismatch between data and model. 



EM algorithm 

The maximization of the likelihood functi on is achieved via th e 
Expectation-Maximization algorithm (EM) dDempster et al. 1977b . 
This algorithm will process iteratively from an initial value of 
the parameters do(b) following a sequence of parameter updates 
8j(b), called 'EM steps'. During the i lh E-step we compute the 
expectation value for the likelihood from the (i - l) th iteration's 
parameters. The i' h M-step maximizes the likelihood (i.e. min- 
imizes the log-likelihood) to compute the i th set of parameters. 
By construction each EM step improves the spectral fit by maxi- 
mizing the likelihood. For a more detail ed review of the spe ctral 
matching EM algorithm used here, see dSnoussi et al. 2001b and 
dDelabrouille et al. 2003b and appendix [A] for the formalism used 
to describe the polarized sky model and data. 

In the MD-MC model presented above there is a scale inde- 
termination on the value of A and R s and only the product AR S A T 
is scale invariant. Thus, to ease the convergence of the algorithm, 
we renormalize each column of A to unity at each EM iteration and 
correct the R s density matrix accordingly so that the product AR S A T 
is unchanged. 



Initialization of the algorithm 

To start-up the EM algorithm the parameters of the fit, 9, need to be 
initialized to reasonable values to avoid exploring local maxima in 
the likelihood function. 

In the case of the mixing matrix, A, we can consider, in a first 
approximation, that the electromagnetic spectrum is the same in 
temperature and polarization for each of the components. There- 
fore, we can concentrate on guessing the electromagnetic spectrum 
in the temperature data where we expect the signal to noise to be 
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larger. When no physical prior is available this can be obtained by 
using the n c dominant eigenvectors of the data density matrix for 
temperature only, R\ T . Each of them represent the change of power 
with frequency for the dominant components in the data. Notice 
that these components and the physical ones are not necessarily 
the same. On one hand, we can have in the data extra components 
which have not been identified as for example residual systematics. 
On the other hand, the electromagnetic spectrum of the physical 
components may present spatial variations as it is, for example, the 
case for the Galactic synchrotron diffuse emission. In the following, 
we will consider that the data contain only identified physical com- 
ponents with spatially constant electromagnetic spectrum. If this is 
not the case, a careful pre-analysis of the initialization parameters 
is needed and this is not discussed in this paper. 

Assuming the mixing matrix previously initialized, the 
physical-components density matrix can be obtained from a noise- 
less fit to the data as follows 

R sB (b) = blockdiag( (A T A o y l A T R y (b)A (A A r )~' ) (7) 
Finally, from A and R s o(b) the noise density matrix is given 

by 

R n ,o(b) = diag( £ ^ (R y (b) - A Q R S . Q A T ) \ (8) 

b 

where n = 2 n b . Here we implicitly assume that the noise is white 
and not correlated between different channels. 



4 SIMULATED OBSERVATIONS 

We have first performed various sets of simulations of the expected 
Planck satellite data to intensively test the algorithm presented 
above. For each of those, we performed 300 realizations consider- 
ing full-sky maps that can contain instrumental noise, CMB, dust, 
synchrotron and free-free. For each realization the CMB and noise 
are changed while dust, free-free and synchrotron are kept un- 
changed. We simulate maps at the LFI and HFI polarized channels, 
30, 40 and 70 GHz for LFI and 100, 143, 21 7 and 353 GHz for 
HFI. These maps are in HEALPix pixelization dGorski et al. 19991) 
and correspond to a 14-month survey. 

(i) [planck a] : We simulate maps at ra sidc = 512 (pixels of area 
~ 47 arcmin 2 ). We include in the simulations CMB emission and 
gravitational lensing. The simulations also contain synchrotron and 
simplified thermal dust emissions as described in section [2] This 
"side permits the reconstruction of the angular power spectra up to 
£ = 1500. The reconstructed spectra will be averaged over bins of 
size 20 in I. With these simulations we test the separation method 
at small angular scales. The simulations contain also instrumental 
noise. 

(ii) [planck b]: We simulate maps at /7 sidc = 128 (pixels of 
area ~ 755 arcmin 2 ) for CMB only. This /7 s j dc permits the recon- 
struction of the power spectra up to ( = 400. This maximum I 
value is enough to study the effect of gravitational waves in the BB 
spectrum which is maximum around ( = 100 for the concordance 
model. The reconstructed spectra are averaged over bins of size 10 
in I. These simulations were performed for fast and intense test of 
the algorithm in the easiest possible case. The simulations contain 
also instrumental noise. 

(iii) [planck c] : We simulate maps at ;i sidc = 128 including 
CMB, simplified thermal dust and synchrotron emissions. The re- 
constructed spectra are averaged over bins of size 10 in I. These 



simulations were performed to check the impact of foregrounds in 
the reconstruction of the CMB power spectra. The simulations con- 
tain also instrumental noise. 

(iv) [planck d] : Same as [planck c] simulations except that 
simplified dust has been replaced by realistic dust and free-free 
emission has been included. 



5 TESTING AND PERFORMANCES ON THE 
SIMPLIFIED-MODEL 

We have applied the MD-MC PolEMICA component separation 
algorithm to the simplified simulations presented above. For each 
set of simulations we have computed the data density matrix R y (b) 
and applied the algorithm with different degrees of freedom: 

(i) First, we assume that the mixing matrix, A, is known. 
We construct the A matrix from the exact electromagnetic spec- 
trum of the anisotropies of each component and fix it in the al- 
gorithm. Therefore, we take as parameters for the fit 6(b) = 
{R s (b), diag(R„(b))} for each bin b. With this test we want to check 
the spatial separability of the components. In the following, we re- 
fer to this type of component separation as A-fixed separation. 

(ii) Secondly, we have performed what we call a CMB semi- 
blind separation. The A matrix is fitted as well as R s and R„. 
but assuming a prior on the CMB electromagnetic spectrum. 
Thus, the columns of the A matrix corresponding to the CMB 
are fixed to unity and the set of parameters for the fit is 8(b) = 
{AjjtcMB,Rs(b), diag(R„(b))} for each bin b. The initialization of A, 
except for the CMB, is performed as described in section [X2l This 
kind of prior in the CMB electromagnetic spectrum is a reasonable 
approximation because the former is well-known. 

(iii) Finally we have performed a blind separation fitting all the 
parameters 6(b) = {A, R s (b), diag(R„(b))} for each bin b including 
the CMB electromagnetic spectrum. 

We have performed these three types of analysis in all the sim- 
plified simulated sets. To ensure the reliability of the results we per- 
formed 10000 EM iterations and checked, for each simulation, the 
convergence of the EM algorithm. In the following we present the 
main results obtained in reverse order going from (iii) to (i) and if 
not stated otherwise we assume white noise R„{b) = o% TEEBB 

5.1 Blind separation analysis 

We present here a blind analysis of the [planck c] simulations for 
which we assume three physical components in the data: CMB and 
simplified-dust and synchrotron emissions. The noise and physical- 
components density matrices as well as the mixing matrix are ini- 
tialized as described in section 13.21 No physical priors are as- 
sumed neither for synchrotron nor dust. For CMB, if not stated 
otherwise, we initialize the electromagnetic spectrum to 1 for tem- 
perature and polarization. This is a reasonable approximation as 
we expect the Planck da ta to be calibrated to better than 1 % 
jPlanck Consortium 20051) . 

Reconstruction of the power spectra 

Figure [T] shows the blind reconstructed CMB temperature and 
polarization power spectra, in y uK 2 , MB , for the 100 GHz channel 
(crosses). The input model is overplotted in blue. The mean and 
the error bars are computed from the analysis of 300 simulations. 
We observe that Cj E , Cj B and Cf B are recovered with no bias, 
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Figure 1. Blind reconstructed CMB power spectra (crosses), Cj T , C EE , C BB , C TE ', C\ B and Cf B , in /JK^MB' at 100 GHz for the [planck c] simulations. 
We overplot in blue the input model. Error bars are computed using a total of 300 simulations. 



up to I = 383, which is the largest accessible C value at n ns ide = 128. 
In the same way, Cf E is accurately reconstructed except for the very 
high ( values for which pixelization problems may appear (see sec- 
tion |5,3l >. We also reconstruct efficiently Cj T although a small bias 
(below 10 %) is introduced at low i mainly due to confusion with 
the synchrotron emission as discussed in the following. The C BB 
spectrum is not recovered at all and a significant bias is observed. 
This bias, due mainly to statistical residual noise as discussed in 
the following section, depends only on the signal to noise ratio and 
does not affect the reconstruction of the other components. 

Figure [2] shows the reconstructed simplified-dust emission 
power spectra, Cj T , Cf E , C BB , Cj E , C] B and Cf B , in pK* MB , at 
353 GHz. For comparison we overplot in red the input model. The 
mean and the error bars displayed were computed using a total of 
300 simulations. The reconstruction is fully efficient for Cj T , Cj E , 
C TB and C EB up to £ = 383. The Cf E and Cf B are accurately recon- 
structed except at I > 350 where a small bias (below 10 %) appears. 
The TB and EB spectra are compatible with zero as expected from 
the input model. These results are consistent with the fact that the 
simulated simplified-dust emission dominate the simulated maps at 
the HFI channels for which the signal to noise ratio is larger. 

Finally, we present in figure [3] the blind reconstructed syn- 
chrotron power spectra, in pK^, MB , at 30 GHz. We overplot in green 
the power spectrum of the input temperat ure and polarization syn- 
chrotron map from dGiardino et al. 20021) . Here again C EE , C BB , 
Cj E , Cj B and Cf B are recovered efficiently. A bias at low I (below 
20 %) is observed for Cj T . This is due, as discussed in the follow- 
ing, to the slight mixing-up of the synchrotron and CMB emissions 
in temperature. 



Reconstruction of the mixing matrix 

The slight mixing up between synchrotron and CMB is better ob- 
served in the reconstructed mixing matrix. The first row of figure 
[4] shows the blind recovered A matrix for dust, synchrotron and 
CMB in the case of the [planck c] simulations. The electromag- 
netic spectrum for T, E and B are respectively traced in green, red 
and blue. For comparison we overplot the input electromagnetic 
spectrum for each of the components (black dashed line). For con- 
venience we have renormalized the electromagnetic spectrum so 
that it is unity at 353 GHz, 30 GHz and 100 GHz for dust, syn- 
chrotron and CMB respectively. It is important to remark that the 
reconstruction of the electromagnetic spectra for temperature and 
polarization is performed independently. We observe that the error 
bars are larger for polarization than for temperature as we would 
expect from the smaller signal to noise ratio in polarization. 

The simplified-dust electromagnetic spectrum is reconstructed 
with no bias both in temperature and polarization even at the low- 
est LFI frequency channels. Furthermore, the reconstructed syn- 
chrotron electromagnetic spectra in polarization are not biased. In 
temperature we observe that the spectrum flattens out at interme- 
diate frequencies between 70 and 217 GHz. This is the cause of 
the slight mixing up between synchrotron and CMB. This mixing 
up does not happen in polarization for which the synchrotron emis- 
sion dominates over the CMB emission. Finally, the reconstruction 
of the CMB electromagnetic spectrum from the T and E modes, al- 
though noisy for the latter, is not biased. However for the B modes 
the reconstruction is very poor because of the very low signal to 
noise ratio (below 10~ 2 for i > 100). 

We have repeated the analysis with no prior in the electro- 
magnetic spectrum of the CMB. We use instead the eigenvector 
corresponding to the third larger eigenvalue of the data density ma- 
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Figure 3. Blind reconstructed synchrotron emission power spectra, Cj T , Cf E , C BB , Cj B , Cj B and Cf B , in /J^mb , at 353 GHz. We overplot in green the 
input model. The error bars are computed using 300 simulations. These are significantly smaller than for the other two components because they do not include 
cosmic variance. Indeed we use the same synchrotron temperature and polarization maps for all the simulations. 
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Figure 4. Reconstructed mixing matrix, A. Columns, from left to right, correspond to the reconstructed A matrix for dust, synchrotron and CMB respectively. 
Rows, from top to bottom, correspond to the blind, blind assuming equal temperature and polarization electromagnetic spectrum and the CMB semi-blind 
analysis of the [planck c] simulations. We display the output electromagnetic spectrum computed for the T (green), E (red) and B (blue) modes as well as 
the input electromagnetic spectrum (dashed black). 



trix. The results for dust and synchrotron remain unchanged. For 
CMB the electromagnetic spectrum at 30 and 353 GHz is not re- 
constructed neither in temperature nor in polarization and the re- 
sults for B are significantly degraded at all frequencies. However, 
the results on the reconstruction of the spatial power spectra re- 
main unchanged for all the physical components including CMB. 
This can be easily understood as the reconstruction of the CMB 
power spectra is mainly dominated by the intermediate frequency 
maps, from 70 to 217 GHz, where the A matrix is accurately recon- 
structed. 



Assuming equal temperature and polarization electromagnetic 
spectrum 

In the previous analysis we have computed the electromagnetic 
spectrum of the physical components independently for each mode 
T, E and B. In a more realistic approach we should consider a sin- 
gle electromagnetic spectrum for the polarization E and B modes 
which may be different from the temperature one. We have re- 
peated the analysis under the above hypothesis and the results re- 
main roughly the same with respect to the reconstruction of the 
spatial power spectra and of the electromagnetic spectrum. For syn- 
chrotron and simplified-dust the E and B modes have roughly the 
same power in our simulations and therefore we expect only varia- 



tions in the error bars. For CMB the E mode largely dominates the 
B mode and therefore we expect no significant contribution from 
the latter to the electromagnetic spectrum reconstruction. 

The differences between the temperature and polarization 
electromagne tic spectra are expected to be small for dust and 
synchrotron i Planck Consortium 20051) and none for the CMB 
dZaldarriaga & Seliak 1997^ Therefore, in the case of a perfectly 
calibrated experiment, we can consider, in a first approximation, 
that the polarization and the temperature electromagnetic spectra 
are the same. The blind analysis of the [planck c] simulations 
under this hypothesis shows no evidence of mixing up between 
synchrotron and CMB. This can be clearly observed in the sec- 
ond row of figure [4] where we represent from left to right the re- 
constructed electromagnetic spectrum of the dust, synchrotron and 
CMB emissions respectively. No bias is observed for any of the 
physical components including synchrotron for which the flatten 
out of the spectrum observed before is not present. 

Figures|5]and[6]show the CMB and synchrotron reconstructed 
spatial power spectra in temperature and polarization. We do not 
observe a bias neither on the synchrotron nor on the CMB TT 
power spectra. Furthermore, we observe that the CMB BB modes, 
although biased at large t values, are fairly reconstructed up to 
f = 60. For the other modes, the results are similar to those pre- 
sented in the previous section. 
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Figure 6. Same as figure[3]but for a blind analysis assuming equal temperature and polarization electromagnetic spectrum. 
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Figure 7. Same as figure[T]but for a CMB semi-blind (black) and A-fixed (red) analysis (see text for details). 



5.2 Semi-blind separation 

The separation method allows us to easily include previous knowl- 
edge on the physical components either as priors or as facts. In 
the previous section we considered a prior on the CMB emission. 
In the following we move a step forward in the analysis assuming 
the CMB electromagnetic spectrum known and performing what 
we call a CMB semi-blind analysis. For this analysis, the columns 
of the mixing matrix corresponding to the CMB are initialized to 
the CMB electromagnetic spectrum and are not updated by the al- 
gorithm. For the other components we consider independent elec- 
tromagnetic spectra for temperature and polarization and they are 
initialized as for the blind analysis. 

The third row of figure [4] shows the reconstructed mixing ma- 
trix for the CMB semi-blind analysis considering the [planck c] 
simulations. The results are similar to those of the blind analy- 
sis discussed before. The dust electromagnetic spectrum is accu- 
rately recovered in temperature and polarization. For synchrotron, 
the polarization spectrum is accurately recovered but the temper- 
ature one flattens out at intermediate frequencies with respect to 
the input model. Therefore, the reconstructed simplified-dust and 
synchrotron spatial power spectra in temperature and polarization 
are similar to those of the blind analysis. The simplified-dust power 
spectra are accurately reconstructed in temperature and polariza- 
tion. For synchrotron the power spectra are also accurately recon- 
structed except for the TT mode which present a slight bias (below 
20%) at large angular scales ({ < 60). In general the error bars are 
smaller for the CMB semi-blind analysis. 

The reconstructed CMB power spectra for the CMB semi- 
blind analysis are represented in black on figure[7] All of them are 
accurately reconstructed with no bias except for the BB mode. For 
the latter the reconstruction is accurate up to £ = 60 and there on 
is biased. This bias is due to residual noise and is not related to the 



uncertainties on the reconstruction of the electromagnetic spectrum 
for the other physical components. To check this we have also per- 
formed a A-fixed analysis assuming the electromagnetic spectrum 
of all physical components known. The results of this analysis are 
overplotted in red on the figure. We observe that reconstruction is 
equivalent to that of the CMB semi-blind analysis but for the error 
bars which are smaller This indicates that the bias in the BB mode 
is mainly due to residual noise as discussed in section \5A\ 



5.3 Reconstruction of the small angular scales 

In the previous section we have fully described the analysis of the 
[planck c] simulations at rc s ;j e = 128 for which the reconstruc- 
tion of the spatial power spectra was limited to i = 383. In some 
cases we have observed small biases in the polarization auto power 
spectra at large ( values which may be due to pixelization prob- 
lems (we exclude in here the bias observed in the CMB BB modes 
which is due to residual noise). To check this hypothesis we have 
also performed the blind, CMB semi-blind and A-fixed analysis on 
the [planck a] simulations for which we can reconstruct the an- 
gular power spectra up to t = 1535. As the resolution of the Planck 
best channels is ~ 5 arcmin a more realistic analysis will require 
simulations at n sjde = 2048 which are far too much time demand- 
ing for our computational capabilities. 

The results obtained for the [planck a] simulations are very 
similar to those for the [planck c] ones. For illustration fig- 
ures [8] |9] and [TO] show the reconstructed temperature and polariza- 
tion power spectra for CMB, simplified-dust and synchrotron in the 
case of a CMB semi-blind analysis. The synchrotron and simpli- 
fied-dust power spectra are accurately reconstructed in temperature 
and polarization. We observe a slight bias in the dust EE and BB 
modes as for the [planck c] simulations but at much larger I val- 
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Figure 8. Same as figure^but for the [planck a] simulations for a CMB semi-blind analysis. 
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Figure 9. Same as figure[2]but for the [planck a] simulations for a CMB semi-blind analysis. 
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Figure 10. Same as figure[3]but for the [planck a] simulations for a CMB semi-blind analysis. 



ues. This indicates that this bias is related to pixelization effects. 
The same effect is observed for the CMB power spectra. The bias in 
the EE modes is present at much larger C than for the [planck c] 
simulations. The reconstruction of the CMB BB modes is accurate 
at low £ (I < 60) and present a residual noise bias at large I values 
as discussed in the previous section. 



5.4 Color noise model 

As seen before, the bias observed in the CMB BB power spectrum 
is most probably due to residual noise from the separation. There- 
fore, it is interesting to check both the accuracy of the noise re- 
construction and the limitations of the white noise model imposed. 
With respect to the latter we have repeated all the analysis presented 
above assuming a color noise model such that the noise power spec- 
tra are estimated for each bin b in £. For the blind analysis the re- 
sults are slightly worse in the sense that the mixing up between 
CMB and synchrotron is more significant. This is not surprising 
since the noise in the data is white and we are artificially reducing 
the number of degrees of freedom in the fit. Actually, the mismatch 
between the data and the physical component power spectra can be 
compensated by changing the noise power spectra. In the case of 
the CMB semi-blind and A-fixed analysis the results for the white 
and color noise model present no significant differences. 

Figure [TT] shows, in black, the reconstructed noise angular 
power spectra, C N t JJ ', Cf ££ , Cf ss , in ^K^, at 100 GHz for the 
[planck a] simulations in the case of a CMB semi-blind analysis. 
We overplot, in black, the power spectra of the noise at 100 GHz 
obtained from 100 realizations of noise-only maps. The noise spec- 
tra are reconstructed down to the accuracy of the estimation of the 
input model both for TT and EE, well below 10~ 3 %. For the 



BB noise spectrum there is a small bias which is of the order of 
5 x 10~ 3 % at t = 200 and around 2 x 10~ 2 % at t = 1500. 

Therefore, to improve the reconstruction of the CMB BB 
modes we need a better estimation of the BB noise power 
spectrum. For this purpose we need to improve the likeli- 
hood maximization algo rithm. For temperature-only separation, 
jDelabrouille et al. 20031) complemented the EM algorithm with 
a direct maximization of the likelihood function via a Newton- 
Raphson algorithm. For polarization similar algorithms can be used 
but due to the degree of complexity of the problem (6 correlated 
modes per physical component instead of 1 in the temperature only 
case) and for the sake of clarity these will be discussed in a forth- 
coming paper. 



6 TOWARDS A MORE REALISTIC SKY MODEL 

After testing intensively our algorithm on our simplified model, de- 
ducting the global performances of the spectral matching recon- 
struction in temperature and polarization we are interested in the 
performances of the algorithm for significantly spatially correlated 
components. For this purpose, we have used the [planck d] sim- 
ulations. 

With this set of simulations, we have performed several types 
of separations. First, we have worked on a joint temperature and 
polarization analysis, similar to the one presented in section[5] Then 
we have considered temperature-only separation, and polarization- 
only separation. For each of the described cases we have applied 
the algorithm with the different degrees of freedom presented in 
section[5] 

Models and recovered data are average over bins of size 10 in 
C beginning at ( = 3. Error bars presented in this section represent 
the dispersion over 100 simulations. Note that rc^e = 128 allows 
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Figure 11. Reconstructed noise angular auto power spectra (black), in fiK^ MB , at 100 GHz for T, E and B for the [planck a] simulations in the case of a 
CMB semi-blind separation. We overplot in red the noise auto power spectra at 100 GHz as obtained from 100 realizations of noise-only maps. 




a theoretical recon struction up to t - 383 in temperature and I ~ 
256 in polarization dGorski et al. 1999b . CMB recovered spectra are 
plotted at 100 GHz while dust at 353 GHz and synchrotron and 
free-free at 30 GHz. 



6.1 Joint temperature and polarization separation 

We have first performed a joint temperature and polarization sepa- 
ration on the realistic model. We present here the results obtained 
in the A-fixed separation case, considering the 4 simulated compo- 
nents. In this particular case, the algorithm is very slow to converge 
and thus 40000 EM iterations have been run. The recovered angular 
power spectra for temperature and polarization are displayed from 



figure [12] to figure Q3] respectively for CMB, realistic dust, syn- 
chrotron and free-free emissions and compared to the input model. 

For the CMB component, on figure[T2]we can see that Cj T and 
Cj E are recovered accurately up to f ~ 400. Recovered spectra for 
Cj B and Cf B are compatible with zero as expected. Cf E is recov- 
ered with a good accuracy up to £ ~ 250 and then is slightly biased 
due to pixelization problems in the HEALPix scheme. Finally, the 
Cf B spectrum is recovered up to t ~ 90 and then is biased with 
residual noise as discussed in section [5741 

The Realistic dust component recovered spectra are displayed 
on figure[l3] We can see that Cj T , Cf E , Cf B , Cj E , Cj B and Cj T are 
recovered with a perfect accuracy up to f = 383. 

The recovered power spectra for the synchrotron component 
are displayed on figure [14] Polarization only power spectra (Cf E , 
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Figure 13. A-fixed reconstruction of the realistic dust C T e T , C EE , C BB , Cj E , Cj B and Cf B power spectra (crosses), in pKj; MB , at 353 GHz in the case of the 
realistic model with 4 components. We overplot in red the input model. Error bars are computed using a total of 100 simulations. 
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Figure 14. A-fixed reconstruction of the synchrotron emission Cj T , Cf E , Cf B , C T t E , Cj B and Cf B power spectra (crosses), in /iK^, MB , at 30 GHz in the case 
of the realistic model with 4 components. We overplot in green the input model. Error bars are computed using a total of 100 simulations. 
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Figure 15. A-fixed reconstruction of the free-free emission C TT , C EE , C BB , C T t E , C TB and C EB power spectra (crosses), in f^^-Q MB , at 30 GHz in the case of 
the realistic model with 4 components. We overplot in orange the input model. Error bars are computed using a total of 100 simulations. 



Cf B and Cf B ) are recovered in good agreement with the input 
model. Cj T is well recovered up to t ~ 200 but then converge 
sharply to a null signal and therefore a residual noise bias similar to 
the one of Cf B occurs. This is not directly visible in the temperature 
spectrum but can be seen in the cross temperature and polarization 
spectra. Indeed, Cj E and Cj B are well recovered up to f ~ 200 and 
further many points are strongly biased. 

Finally, results corresponding to the free-free emission com- 
ponent are displayed in figure [15] Cj T is recovered with a good 
accuracy up to I = 383. As for polarization the input signal is null, 
we can see in the recovered spectra the overall behavior of our al- 
gorithm described in section|5] 

The algorithm has also been run for CMB-fixed and Blind sep- 
arations. For both of them, excepting the dust component which 
is well constrained in all cases due to its dominant power at high 
frequencies, the algorithm fails to converge and then components 
are mixed and results strongly biased. This may come from the fact 
that free-free and synchrotron electromagnetic spectra are similar 
and to the fact that all the Galactic emissions have strong spatial 
correlations. In the following section we will address this problem 
and show that this mixing that prevents the convergence of the algo- 
rithm is mainly due to the separability problem that occurs in tem- 
perature. In section [631 we will see that the separation performed 
on sets of Q and U maps has not this separability problem. 

Notice that when not considering the free-free emission in the 
simulations, results are very similar and the same performances of 
the algorithm with respect to the level of prior we assume are ob- 
served. 



6.2 Temperature-only separation 

We have performed a temperature-only separation on the realis- 
tic model. For this we consider sets of / maps and the algorithm 
solve the spectral matching equations for T modes only, like in the 
SMICA algorithm jDelabrouille et al. 20031) . We present here the 
results obtained in the A-fixed separation case, considering the 4 
simulated components. 

Recovered Cj T spectra are displayed in figure [T6] for syn- 
chrotron, realistic dust, CMB and free-free. We can see that except 
the synchrotron spectrum which start to be biased at I ~ 200, for 
the reason advanced in the last section, spectra are recovered with 
a good accuracy for the dust, CMB and free-free components up to 
t = 383. We have also performed the separation for the CMB-fixed 
and the Blind cases. As for the joint analysis, excepting again the 
dust component for which the spectrum is recovered efficiently, the 
algorithm fails to converge. Note that in the literature, no method 
has successfully separated synchrotron from dust on a noisy simu- 
lated mixture of CMB and astrophysical foregrounds, working on 
all sky maps. 

Comparable results and the same performances of the algo- 
rithm with respect to the level of prior we assume are observed 
when no free-free component is included in simulations. 



6.3 Polarization-only separation 

Finally, we have performed a polarization-only separation on the 
realistic model. For this we consider sets of Q and U maps and 
the algorithm solve the spectral matching equations for E and B 
modes (allowing reconstruction of Cf E , Cf B and Cf B ). Notice that 
as we suppose that the free-free emission is not polarized, the sets 
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Figure 16. A-fixed temperature-only reconstruction of the Cj T power spectra (crosses) in fiK^ MB in the case of the realistic model with 4 components. Results 
are displayed from left to right for the 4 different components, i.e synchrotron (at 30 GHz), realistic dust (at 353 GHz), CMB (at 100 GHz) and free-free (at 
30 GHz). Input models are overplotted in solid lines. Error bars are computed using a total of 100 simulations. 
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Figure 17. A-fixed polarization-only reconstruction of the CMB Cf E , Cf B and Cf B power spectra (crosses) in pK^. MB at 100 GHz in the case of the realistic 
model with 3 polarized components. Input models are overplotted in blue. Error bars are computed using a total of 100 simulations. 



of polarized maps used here only contain CMB, realistic dust, syn- 
chrotron and noise. 



case, only the error bars are larger. However, the Cf spectrum is 
not recovered. 



CMB semi-blind separation 

We have performed a CMB semi-blind separation on the 
polarization-only set of maps. In this case, the CMB electromag- 
netic spectrum is initialized to 1 and kept fixed (see section|5j. 

The results of this analysis on the CMB power spectra are dis- 
played in figure [T7] We can see that Cf E are reconstructed without 
bias up to t ~ 300, as expected with this pixelization scheme and 
that Cf B is compatible with zero. Cf B is efficiently recovered up to 
I ~ 50 and then is biased with residual noise from the separation. 

The realistic dust and synchrotron emissions recovered power 
spectra are displayed in figure [T8] and [T9] In both cases we can 
see that Cf E , C BB and C EB are reconstructed very efficiently up to 
£ = 383. 



Blind separation 

We have also performed, for the polarization-only separation, a 
Blind separation. The results of this analysis for dust and syn- 
chrotron are not displayed, but for both, Cf E , C BB and Cf B are re- 
constructed very efficiently up to £ = 383. The reconstruction of 
the CMB component is displayed in figure [20] For C EE we can see 
that the reconstruction is similar to the one in the CMB semi-blind 



6.4 Discussions 

From the previous analyzes we have clearly identified a separa- 
bility problem when dealing with more than one realistic diffuse 
Galactic emission component. This problem appears both in the 
joint temperature and polarization and in the temperature-only ana- 
lyzes, but not in the polarization-only one. This would indicate that 
it is mainly due to the high level of correlation of the Galactic dif- 
fuse emission in temperature both in the Galactic plane and at high 
Galactic latitudes. The current version of our algorithm assume un- 
correlated components and therefore we expect it to behave badly 
when they are correlated. Work is in progress to adapt PolEMICA 
to account for spatially correlated components. 

For polarization-only separation the correlation problem 
seems to be not significant and the CMB semi-blind and blind an- 
alyzes are possible. Therefore, for this case, we can evaluate the 
loss of accuracy in the reconstruction of the CMB signal due to the 
foreground contamination. For this, we compare the error bars on 
the CMB power spectra for the A-fixed, CMB semi-blind and blind 
separations to the one obtained in the case of a A-fixed separation 
on [planck b] simulations containing CMB and noise only. 

Results are presented in figure [2T1 C EE , Cf B and C EB recon- 
struction error bars behave similarly with respect to the different 
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Figure 18. A-fixed polarization-only reconstruction of the realistic dust Cf E , C BB and C EB power spectra (crosses) in juK^g at 353 GHz in the case of the 
realistic model with 3 polarized components. Input models are overplotted in red. Error bars are computed using a total of 100 simulations. 
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Figure 19. A-fixed polarization-only reconstruction of synchrotron emission C EE , C BB and Cf B power spectra (crosses) in /iKj; MB at 30 GHz in the case of 
the realistic model with 3 polarized components. Input models are overplotted in green. Error bars are computed using a total of 100 simulations. 



algorithm priors. For each of them, A-fixed and CMB semi-blind 
error bars are of the same amplitude and are between 25 and 50 % 
larger (respectively a factor 1.26 for Cf E , 1.50 for Cf B and 1.29 for 
Cf B ) than the reference A-fixed CMB only case. This means that in 
the context of our realistic model, we have no need to put priors on 
the foregrounds electromagnetic spectra in polarization to perform 
an efficient separation. In the Blind case, error bars of the recon- 
struction are increased by a factor 1.59 for Cf E , 13.5 for C BB (but 
the reconstruction is biased) and 2.52 for Cf B . On the other hand, 
these results stress the fact that having priors on the CMB electro- 
magnetic spectra inside of our algorithm helps to perform a more 
accurate separation. 



7 SUMMARY AND CONCLUSIONS 

We present in this paper the PolEMICA algorithm which 
is an extension to polarization of the SMICA temperature 
MD-MC blind compon ent separation method developed by 
dDelabrouille et al. 20031) . Both algorithms work in harmonic space 
and are based on the spectral matching of the data to a noisy linear 
mixture of uncorrected physical components using the EM algo- 
rithm to maximize the likelihood function. By contrast to the tem- 
perature data which are described by a single scalar quantity T, the 
combined temperature and polarization data are described in har- 
monic space by three correlated scalar quantities T, E and B cor- 
responding to the /, Q and U Stokes parameters in real space. We 
have developed a new formalism to jointly deal with the 6 resulting 



auto and cross angular power spectra, TT, EE, BB, TE, TB and 
EB. 

Using this formalism we have constructed the likelihood 
function and proved that the EM algorithm can be also applied to 
polarization data. Under the assumption of uncorrelated Gaussian 
distributed components and noise, the free parameters in the fit 
are the mixing matrix describing the electromagnetic spectrum 
of the physical components for T, E and B, the temperature and 
polarization angular power spectra of the physical components and 
the temperature and polarization noise power spectra for each of 
the detectors. 

We have, as a first approach, intensively and successfully 
tested the PolEMICA method on simulations of the Planck satel- 
lite experiment considering a 14-months nominal mission and no 
systematic effects. For these tests, we suppose a simplified lin- 
ear model for the sky emission including CMB, synchrotron with 
constant spectral index and simplified-dust (Gaussian realization) 
emissions. We construct full sky maps for all the polarized chan- 
nels from 30 to 353 GHz including at least one of the above physi- 
cal components and considering white noise and infinite resolution. 

The method permits blind separation on these simulations 
allowing us to reconstruct the noise and physical component's 
temperature and polarization power spectra as well as the mixing 
matrix when we consider equal electromagnetic spectrum in T, E 
and B. When we relax this hypothesis the reconstruction of the 
electromagnetic spectrum for the CMB B modes is significantly 
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Figure 20. Blind polarization-only reconstruction of the CMB C EE , C BB and Cf B power spectra (crosses) in /jK^. mb at 100 GHz in the case of the realistic 
model with 3 polarized components. Input models are overplotted in blue. Error bars are computed using a total of 100 simulations. 
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Figure 21. Error bars of the reconstruction of the CMB C EE , C BB and C EB power spectra in /iK^, MB at 100 GHz in a Blind (blue), CMB semi-blind (green) 
and A-fixed separations of a 3 components mixture containing CMB, realistic dust and synchrotron and in a A-fixed separation of a set CMB and noise only 
simulations (black). 



degraded as could be expected because of the low signal to noise 
ratio. These results indicate that the PolEMICA method allows 
us to both constrain the electromagnetic spectrum of the physical 
components and also to inter-calibrate the data based on the 
reconstructed CMB electromagnetic spectrum. 

After setting the general performances of the algorithm, we 
have performed the separation on a more realistic model that in- 
cludes realistic-dust, synchrotron and free-free components in sec- 
tion^ We have encountered in this case a separability problem, that 
mixes up components and prevents the algorithm to converge, when 
performing blind separations. We have shown that this is due to 
spatial correlations between Galactic components in temperature. 
Thus, when working on sets of Q and U maps and maximizing the 
likelihood for EE, BB and EB modes only, this separability prob- 
lem does not occur and CMB semi-blind and blind separations are 
possible. For this polarization-only case, we have shown that con- 
sidering our realistic sky model and our algorithm, in the Planck 
case, we have no need to put priors on the Galactic components 
electromagnetic spectra to reconstruct the CMB polarized power 
spectra. Nevertheless, adding priors on the CMB electromagnetic 
spectrum helps to perform a more accurate separation. 

Finally, real experiments present finite resolution, partial ef- 
fective sky coverage, systematic effects and, often, correlated noise. 
All these issues must be dealt with by the component separation al- 
gorithms and will with no doubt significantly limit the precision 



to which the CMB signal may be reconstructed. PolEMICA, as it 
was already the case for SMICA, can account for beam and fil- 
tering smoothing. Systematic effects and correlated noise can be 
modeled as extra components in the data for which the spectral 
dependence can be estimated in a blind analysis. Moreover, the 
strong spatial correlation in temperature between Galactic physical 
emissions: dust, synchrotron and free-free, is a major problem for 
blind component separation algorithms which generally assumed 
uncorrelated components. Although not observed yet, we can also 
imagine spatial correlation of the Galactic emissions in polariza- 
tion. Work is in progress to adapt the PolEMICA algorithm to the 
case of correlated components. 

In addition, foreground emissions have in general spatially 
varying electromagnetic spectra far beyond the simple linear model 
presented here. Work is also in progress to adapt PolEMICA to the 
case of foregrounds with spatially varying electromagnetic spec- 
trum. 
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APPENDIX A: MD-MC POLARIZATION SKY MODEL 

We discuss here the formalism developed to describe the temperature and polarization observations as a noisy mixture of independent 
components. 

In the following we assume full sky observations at two frequencies v\ and v 2 and a simple linear model for the sky emission with two 

components ci and c 2 . 

In this case, equation l[2j reads 
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where y '^ and sT 1 for X = {T,E,B) are the coefficients of the spherical harmonic decomposition of the input sky observations and of the 
components of the sky model respectively. The coefficients f* c , correspond to the electromagnetic spectrum of the component cj at the 
frequency of observation v, . Note that the mixing matrix, A, has dimensions 3 • riy X 3 • n c . 



We define the noise, R„(b), signal, R s (b), and data, R y (b), density matrices used in equation l(5j for each bin, b, as follows 



R : (b) = — Y< Z(t)z(lf > , Z = {«, s, y) 

hx> b 

Assuming the noise uncorrelated between detectors the noise density matrix is diagonal 
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In the same way assuming independent physical components the signal density matrix is block diagonal and reads 
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Finally the data density matrix can written by blocks as follows 
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